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A SIMPLIFIED FORMULA FOR IONOSPHERIC FARADAY 
ROTATION AT FREQUENCIES ABOVE 100 MC/S 


T. J. Lynch 

Advanced Development Division 


ABSTRACT 

A simple formula for ionospheric Faraday rotation 
is derived starting from Maxwell's equations and using 
the following major assumptions: the ionosphere is a 

homogeneous gyrotropic medium with a scalar permea- 
bility and a tensor permittivity; the earth's magnetic 
field is constant and directed along the path of propaga- 
tion; the ionosphere imposes no attenuation on the linearly 
polarized entering wave; and the carrier frequency is of 
the order of 100 Mc/s or higher. 
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A SIMPLIFIED FORMULA FOR IONOSPHERIC FARADAY 
ROTATION AT FREQUENCIES ABOVE 100 MC/S 


INTRODUCTION 

When only an approximate estimation of ionospheric Faraday rotation is re- 
quired for frequencies above 100 Mc/s, it is possible to derive a simplified 
formula for this purpose by making a number of assumptions. These assump- 
tions lead to a model of the ionosphere which is idealized, and which is not an 
absolute representation of the spatial and temporal characteristics of the actual 
ionosphere. However, this model does permit rapid computation of an approxi- 
mate value of Faraday rotation as a function of frequency, average ionospheric 
electron density, and average ionospheric magnetic flux density. 

This paper presents the related assumptions, derivation, and use of this 
simplified formula. 


ASSUMPTIONS 

The ionosphere is assumed to be a gyrotropic medium. A gyrotropic medium 
is defined as an anisotropic medium in which the permittivity tensor is not sym- 
metric. In the present case, the finite earth's magnetic field acting on the 
ionospheric charged particles brings about this asymmetry in the permittivity 
tensor. 

The earth's magnetic field vector is assumed to be directed along the prop- 
agation path, and to have a constant magnitude along this path. 

The ionospheric electron density is assumed to have a constant value along 
the propagation path. 

The radiated wave entering the ionosphere is assumed to be a linearly 
polarized TEM (transverse electromagnetic) wave. 

The carrier frequency is assumed to be in the order of 100 Mc/s or 
higher. 

The permeability of the ionosphere is assumed to be a scalar quantity equal 
to the permeability of free space. 
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The ionosphere is assumed to exert no attenuation on the propagating wave. 


SIMPLIFIED FORMULA 

By use of the above assumptions , the following simplified formula for one- 
way ionospheric Faraday rotation may be derived as shown in the Appendix. 
Approximately, the per-unit-path-length Faraday rotation is given by: 

2 

BN 

10 3 in radians/km (1) 

z 2 e 0 co? 

where in the rationalized meter-kilogram-second (MKS) unit system: 

e is electronic charge = 1.6 x 10~ 19 coulomb 

m is electronic mass = 9.1 x 10” 31 kilogram 

e Q is permittivity of free space = 8.85 x 10" 12 farad per meter 

c is velocity of light = 3 x 10 8 meters per second 

co is carrier frequency in radians/sec 

B is magnetic flux density in webers/m 2 

N is electron density in electrons/m 3 

Combining constants, the working formula in terms of f (cycles per sec) 
becomes: 



n 

z 


1.355 x 10 9 



in degrees /km 


(la) 


RANGES OF VALUES FOR B AND N 

The assumptions made in developing (1) involve a constant value for N and 
a constant value of B along the propagation path. In addition, the direction of 
the earth's magnetic field is assumed to be identical to the direction of 
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propagation. These rather unrealistic assumptions require that the choices of 
the values of N and B be made with these assumptions in mind. 

The value of B to be used in (1) should be an average one for that region of 
the earth’s surface involved. World maps 1 of the earth’s magnetic field intensity 
and direction can be used to pick a value of B for estimation of Faraday rotation. 
It is significant to note that the intensity of the magnetic field does not vary 
widely over the earth, and that an average value of 0.5 x 10" 4 webers/m 2 is a 
good starting point in the absence of more detailed information, as shown in the 
footnote in the Appendix. Since it is actually the component of the earth's mag- 
netic field colinear with the direction of propagation that gives rise 2 to the Fara- 
day rotation given by (1), the above assumption of a completely colinear magnetic 
field will result generally in an upper-bound type of estimate of Faraday rotation. 

The value of N to be used in (1) depends upon the thickness of the homo- 
geneous spherical shell one wants to assume for the model of the ionosphere. 

It also depends upon the point in the sunspot cycle, and whether the local condi- 
tion is day or night. As indicated in the footnote in the Appendix, for a shell 
that is assumed to range from 100 km to 600 km above the spherical earth, the 
daytime value of average electron density N varies 1 from 0.4 x 10 12 to 10 12 
electrons/m 3 as the 11-year sunspot cycle goes from minimum to maximum. 

For a thinner shell model, correspondingly larger values of N must be used to 
be consistent with a constant value of the integral J Ndh , where h is height above 
the earth's surface. 
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APPENDIX 


DERIVATION OF THE SIMPLIFIED FORMULA 


A linearly polarized TEM wave entering a gyrotropic medium will undergo 
a rotation of field components, and it is this rotation, as a function of the 
parameters of the medium and the frequency, that we wish to formulate. 

A useful starting point 3 is the representation of the linearly polarized wave 
by the sum of two circularly polarized waves rotating in opposite directions. If 
a right-handed coordinate system of x, y, and z, with propagation in the positive 
z direction is used, then at z = 0, the electric field vector is 


E = a E. 

X 


( 2 ) 


Let a counter-clockwise-rotating circularly polarized wave be given by: 


v - E - • E 

E = a — + a i — . 

ccw x 2 y 2 


(3) 


Let a corresponding clockwise-rotating circularly polarized wave be given 


by: 


^ _ E _ . E 

i, = a — - a j — 

cw x 2 y 2 


(4) 


Then 


E = E + E . 

CCW CW 


(5) 
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The Faraday rotation will be found from the phase changes in these two ro- 
tating waves as they propagate through the gyrotropic medium. 

Applying Maxwells' equations to each rotating field separately, for the 
counter-clockwise wave 


BH 


VXE 


= -Mo 


Bt 


( 6 ) 


Taking the curl of both sides, 


V X (V X E 


B(V x H ) 

. ccw 7 

») = - Mq dt 


(?) 


Now instead of using 


V X H = i + e n 

ccw 0 


BE 

CC 

~Bt" 


( 8 ) 


use 


BE 

V X H O) 

CC w Bt 


where e is a permittivity tensor which takes into account the interaction of the 
magnetic field of the medium with the charged particles therein. 

But 


VX(VXEccw) =?(V -E ccw ) - V 2 E c cw’ 


( 10 ) 


6 



where 


V 


E 

ccw 


3(E ) 3(E) 3 (E ) 

v ccw 7 x ^ v ccw 7 y v ccw / z 

3 x 3 y 3 z 


(ID 


Now since E ccw is a field component of a TEM wave, there is no z component, 
and there is no spatial variation of E in the x or y direction. Thus V • E„„,„ 

CCW CCW 

= 0. Then 


V 2 E 


3 2 E 


= f*o 


3t 2 


( 12 ) 


Assume a sinusoidal time variation for E in (3); 


E = E 0 ei" 1 . 


(13) 


Substituting (13) and (3) into (12), 


V 2 E = -a> 2 

ccw 


i E 


(14) 


Assuming the external magnetic field to be oriented in the positive z direc- 
tion, the permittivity is given by: 


e = e„ 


€ 11 € 12 


£ 21 € 22 


0 


'33 


(15) 
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with 



C 


( 16 ) 



C 



(18) 


where 


a) p is the plasma frequency = y^Ne/e Q ) (e/m) 
co is the angular cyclotron frequency = (e/m) B. 

c 

When the counter-clockwise wave propagates through the gyrotropic medium, 
the expression for each field vector takes on a factor e~^ccw z . Since there is no 
spatial variation in the x and y directions, 


V 2 E 


ccw 


21 

Bz 2 


(E ). 

' ccw 


(19) 


Then (14) may be written as three equations: 


21 
Bz 2 


E/2 


£ U e i2 0 


E/2 

j E/2 

“7 z 0 

e ccw = - co 2 e 0 

£ 21 e 22 0 


j E/2 

0 


_ 0 0 e 33_ 


0 


■y z 

' CCW 


( 20 ) 
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Differentiating either the first or second equation in (20) and dividing through 
by E/2 e ^ccw 


(/cc.) 2 =- co f^o e o ( £ n + i e i2^ 


(21) 


Assuming no attenuation in the wave, then y = j f3, and 


^ccw = + j £ 12 )1/2 - 


( 22 ) 


By parallel procedure for the clockwise wave: 


^ " 1^0 € 0 ( £ ll - j € 12 )1/2 - 


(23) 


A, • 


Now it remains to express the Faraday rotation angle in terms of /3 and 


From (5) for any value z that places the wave within the gyrotropic medium, 
the electric field vector can be written: 


-j/3 Z — “ j f3 Z 

+ E 

c cw cw 


E = E e ,rcc *' + e 


E = a E + a E . 
xx y y 


(24) 

(25) 


From (3) and (4) 


E = — (e 
* 2 


E , -i£ CC w z , 


+ e ) 


l? E ~i ^ccw^ • _ j^cw Z N 
E r_ (j e cc * - J e ) 

y 2 


(26) 

(27) 
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Now the Faraday rotation angle is given by 


n = tan -1 — . 

E 

X 


( 28 ) 


Factoring out e~ ,/2 (/3ccw + /3cw>Z from (26) and (27) obtains 


E = Ee 


-j/2 (/3 CCW + /W Z 


f P -J3 

* r*u/ 1 i 


cos 


cw ' ccw 


(29) 


- j /2 (8 + 8 )Z I P nw " 

E - -Ee J ^ccw ^cw ; 


cw ' ccw , 

sin l r 1 ^ 


(30) 


then 


P - /3 , 

- sin l ) z 

"y 

E //3 - /f 

x I r cw ' ccw , 

COS ( - 1 Z 


or 


then 


5 fP - y8 \ 

y .= - tan (— 1- — - c — | z 


E 


(31) 


n = P ccw ) z . 


(32) 


Substituting (22) and (23) in (32): 


fi GJ / 

7=2 ^ 


e 0 [(e il + j e i2^ 1/2 - (e il ' j £ 12' )1/2] 
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or 


n 

z 


CO / 

— =7T / M n 6 


1 1 


1/2 


1 + j 


12 


'll, 


- 1 


e \ 1/2 ‘ 


(33) 


From (16) and (17) 


12 


ii 


0^ 

p 


(i) 


co 2 - a? \ ex? 

c p 


or 


12 


ii 




/cd \ 2 r 

\C*> I \ Ci 

\ 9 / \ 


"’ A 2 


- 1 


(34) 


If (34) is computed for a frequency of 100 Mc/s and some "typical" values* of 
N = 1 0 12 and 3 = 0.5 * 10 -4 , then = 8.8 * 10 6 ; cu = 5.66 x JO 7 ; co = 6.28 x 10 8 
and j (e 12 /e n ) = 1.15 x 10" 4 . 


*From Johnson * , the value of B ranges from 0.25 x 10 4 to 0.70 x 10~ 4 on the total intensity 
plot on page 121. Also from Johnson 1 , on page 28, the integral of the electron density height 
profiles for daytime can be approximated as follows: 

For minimum of sunspot cycle: I Ndh £ 2 x 10 14 - - ^ . C ?/ on S - km 


For maximum of sunspot cycle: | Ndh £ 5 x 10 14 e * ec *r° ns _ 

J m 3 


Then for a 500-km thick ionospheric shell with uniform electron density from 100 to 600 km above 
the earth: 

For minimum of sunspot cycle: N £ 0.4 x 10 12 e ^ ect ^ ons 

m 3 


For maximum of sunspot cycle: N £ 10 


12 electrons 
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The binomial expressions in (33) may then be expanded in series form, retaining 
just the first two terms of each series to a good approximation, as follows: 



1 ; €l2 


(35) 


Then 



n 


'"Vv 

'■v 


z 


^0 e o g ll 
2 



(36) 


(37) 


Substituting for e 1 1 and e 12 , 


_ L 

Z 2 [(* 1 - O) 2 + a,J) (0)2 -c 2 )] 1/2 


(38) 


Note that the above approximation improves as the frequency increases. 
Now at the "typical" values for N and B discussed above, compare o> 2 , a; 2 
and co 2 for a frequency of 100 Mc/s. 


CO 2 a 7.74 X 10 13 ; co 2 = 3.20 x 10 ls ; o 2 . 3 94 x 10 i7 

C P 


Therefore for these values: 


2 2 2 

00 z << 00 z << 00 z . 

c p 
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Then (38) becomes 




0 ( 39 ) 
— ^ — • 

z 2co 2 


Substituting for a> and a> c and letting V ji Q e Q = l/c, the simplified formula be- 
comes 


a _ e (|) 2 BN 

z — t 10 3 radians/km. (40) 

z 2 e 0 c oj 2 
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